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ABSTRACT 

We present a combined X-ray and optical analysis of three bimodal galaxy clusters selected as merging candidates at z ~ 0.1. These 
targets are part of MUSIC (MUlti-Wavelength Sample of Interacting Clusters), which is a general project designed to study the 
physics of merging clusters by means of multi- wavelength observations. Observations include spectro-imaging with XMM-Newton 
EPIC camera, multi-object spectroscopy (260 new redshifts), and wide-field imaging at the ESO 3.6m and 2.2m telescopes. We build 
a global picture of these clusters using X-ray luminosity and temperature maps together with galaxy density and velocity distributions. 
Idealized numerical simulations were used to constrain the merging scenario for each system. We show that A2933 is very likely an 
equal-mass advanced pre-merger ~ 200 Myr before the core collapse, while A2440 and A2384 are post-merger systems (~ 450 Myr 
and ~ 1.5 Gyr after core collapse, respectively). In the case of A2384, we detect a spectacular filament of galaxies and gas spreading 
over more than l/r ! Mpc, which we infer to have been stripped during the previous collision. The analysis of the MUSIC sample 
allows us to outline some general properties of merging clusters: a strong luminosity segregation of galaxies in recent post-mergers; 
the existence of preferential axes -corresponding to the merging directions- along which the BCGs and structures on various scales 
are aligned; the concomitance, in most major merger cases, of secondary merging or accretion events, with groups infalling onto the 
main cluster, and in some cases the evidence of previous merging episodes in one of the main components. These results are in good 
agreement with the hierarchical scenario of structure formation, in which clusters are expected to form by successive merging events, 
and matter is accreted along large-scale filaments. 

Key words. Galaxies: clusters: individual: Abell 2384 - Abell 2440 - Abell 2933 - X-rays: galaxies: clusters - galaxies: intergalactic 
medium 



1. Introduction 

In standard Cold Dark Matter cosmological models, including 
the concordance ACDM, the general growth of structures starts 
from the primordial density fluctuations generated by inflation, 
and is driven by gravity in a hierarchical way, i.e. smaller struc- 
tures form first, then merge into progressively more massive 
systems; however, if the acceleration of the expansion is due 
to a cosmological constant, the process of structure formation 
will completely stop in the future (Krauss & Starkman 2000; 
Nagamine & Loeb 2003; Busha et al. 2005). 

Galaxy clusters are the largest gravitationally bound objects 
of the hierarchy; they accrete smaller groups coming from the 
filaments of the cosmic web and occasionally merge with other 
clusters of comparable mass, releasing an exceptionally high 
amount of energy. Merging clusters are therefore ideal labora- 
tories to study the process of structure formation and how its 
affects galaxy evolution. In this scenario, one expects to find 
a large fraction of clusters in their formation process at high 
redshift (which is corroborated by the large fraction of irregu- 
lar morphologies observed). However, observations of clusters 
at high redshift are difficult and time consuming, so an alter- 
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native choice is to search for rarer but more easily observable 
merging candidates at low redshift. In this way, we can probe in 
detail the merging signatures and try to shed some light on the 
process of cluster formation and how it affects galaxy evolution. 

A crucial progress in the study of merging clusters was pro- 
vided by the spectral and imaging capabilities of the last gen- 
eration of X-ray satellites. Before then, X-ray spectroscopic in- 
formation had not been available or of very poor quality, and the 
main information about mergers was based on morphology. With 
the spatially resolved spectroscopy and high resolution imag- 
ing offered by Chandra and XMM, the situation has radically 
changed, and sophisticated algorithms have been developed to 
achieve good and reliable temperature maps (see for instance 
Bourdin et al. 2004), since the temperature is the more accu- 
rate tracer of the energy transfer from the collision to the X-ray 
gas itself. Strong signatures of the merging events have been de- 
tected in these maps (Vikhlinin et al. 2000, Markevitch et al. 
2000 and 2002), and cold fronts and bow shocks are now well 
established as common merger features. Thanks to this observa- 
tional progress, major mergers of galaxy clusters now appear far 
more complex than previously foreseen. 

The full understanding of the complex processes at work 
in merging requires dedicated numerical simulations. Much 
progress has been made in this field, starting from the pioneer- 
ing works of Schindler and Bohringer (1993) and Roettiger et 
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al. (1997). Ricker and Sarazin (2001) (hereafter RS01) described 
the violent relaxation of gas in a dark matter potential well and 
a variety of idealized merging systems, paying special attention 
to the impact parameter and the mass ratio between units. Poole 
et al. (2006) analyzed merging of idealized relaxed clusters with 
sophisticated SPH simulations including cooling and star for- 
mation, and detected the major transient signatures existing in 
observed temperature maps. 

Combining optical with X-ray data has been shown to be 
extremely effective in unveiling the complex history of merging 
clusters (to mention a few, Arnaud et al. 2000 and Maurogordato 
et al. 2000; Donnelly et al. 2001; Barrena et al. 2002; 2007; 
Boschin et al. 2004; Owers et al. 2009). These studies have re- 
vealed various peculiar properties of the galaxy distribution in 
the individual merging clusters, such as strong signatures in the 
density and velocity distribution, and strong alignments effects. 
However, a larger sample with both X-ray and optical observa- 
tions is clearly needed to test the generality of these properties, 
and their dependence on the merging stage. 

Motivated by these reasons, we started an observational 
program, MUSIC (MUltiwavelength Sample of Interacting 
Clusters) to define and analyze a sample of merging clusters at 
different stages of the merging process. The first targets (A3921, 
A1750, A2065, and A2255) were selected from the merging 
clusters observed during the XMM Guaranteed Time by one of 
us (J.L. Sauvageot), and we established an optical follow-up pro- 
gram at ESO. 

To cover systems representing various stages of merg- 
ing, we extended the initial sample by including targets from 
Kolokotronis et al. (2001), and compared the density distribu- 
tions of the gas (from ROSAT/HRI maps) and of the galaxies 
(from APM maps). We selected the systems with a clear bimodal 
morphology in both the X-ray and the optical maps as proba- 
ble pre-merger candidates, and those with a distorted morphol- 
ogy and a pronounced segregation between gas (collisional) and 
galaxies (non collisional), i.e. with features that are expected of 
mergers at a more advanced stage. We selected clusters around 
z ~ 0.1, because at this redshift evolutionary effects on galaxies 
are negligible and the field of view of X-ray and optical instru- 
ments offers a good coverage of the system (30' with XMM and 
wide field imagers such as WFI at the ESO 2.2m, corresponding 
to 2.3 h" 1 Mpc at z = 0.1). With a 4m class telescope, spec- 
troscopy is also easily performed as faint as R < 19, a limit 
corresponding to ~ L* + 3.5 at this redshift (with a best- fit pa- 
rameter M* = -22.97 estimated by Popesso et al. (2005), fitting 
a Schechter function to the cluster galaxy luminosity function of 
the RASS-SDSS cluster survey). The aims of the program are to 
characterize the merging scenario (epoch of collision, geometry, 
and mass ratios) and to test the impact of the merging process on 
the properties of galaxies, in particular the star formation. For 
example, the combined X-ray and optical analysis allowed us to 
determine the merging scenario of the first two clusters (a pre- 
merger for A 1750, Belsole et al. 2004, and a recent post-merger 
for A3921, Belsole et al. 2005, Ferrari et al. 2005). 

Here we present the analysis of three additional merging 
candidates: A2933, A2440, and A2384, based on XMM X-ray 
spectro-imaging obtained in 2005 and 2006 (A2933 and A2440), 
or retrieved from the XSA database (A2384) and optical obser- 
vations carried out in 2003, 2005, and 2007 at ESO, including 
wide-field multi-band imaging with the Wide Field Imager on 
the 2.2m telescope, and multi-object spectroscopy with EFOSC2 
at the 3.6m telescope. 

The clusters A2933, A2440 and A2384 share common fea- 
tures: an irregular morphology (they have been classified as 



Bautz-Morgan Type III, II, and II/III, respectively), low Abell 
richness classes (1, 0, and 1), and redshifts slightly less than 
0.1 (0.0925, 0.0906, and 0.0943 according to Struble and Rood 
1999). In the optical, A2933 was observed as a target of the ESO 
Nearby Abell Cluster Survey (Katgert et al. 1996) and the south- 
ern Abell redshift survey (Muriel et al. 2002; Way et al. 2005; 
Coenda et al. 2006). A2440 was identified as a pre-merging clus- 
ter from the dynamical and X-ray/optical analysis of Beers et al. 
(1991) and Mohr et al. (1996). A2384 is also a classical exam- 
ple of a bimodal cluster (Ulmer and Cruddace 1982; West, Jones 
and Forman 1995). A weak-lensing analysis of A2384 was per- 
formed by Cypriano et al. (2004). Both A2440 and A2384 be- 
long to the flux-limited sample of bright clusters of galaxies from 
the southern part of the ROSAT All-Sky Survey (de Grandi et al 
1999). 

In Sects. 2 and 3, we present the data, the reduction pro- 
cedure, and the methodology used in X-ray and optical respec- 
tively. In Sects. 4, 5, and 6, we analyze the data for A2933, 
A2440, and A2384 respectively and propose a merging sce- 
nario for each cluster. In Sect. 7 we present our numerical 
simulations of merging clusters that are used in Sect. 8 to re- 
fine the scenarios. In the following, we adopt the cosmologi- 
cal parameters of a ACDM model with Q M =0.3, £2a=0.7, and 
Hq = 70 km s" 1 Mpc -1 . With these parameters, at z = 0.1 one 
degree corresponds to a physical length of 6.6 Mpc. 



2. X-ray observations: scientific products 

The three systems presented in this paper were observed by 
XMM-Newton EPIC cameras for about 30 ks in full frame mode 
with the medium filter. The observation were screened for pro- 
ton flares on a high and low energy light curve basis. A summary 
of the effective exposure time remaining after this screening pro- 
cess is provided in Table [T] 

The brightness of the X-ray emitting intra-cluster medium 
(ICM) was mapped from a multi- scale algorithm using Haar 
wavelets, especially suited to denoise images dominated by shot 
noise. First proposed by Jammal & Bijaoui (2004), this algo- 
rithm was adapted to analyze X-ray astronomical images and 
correct signal distorsions related to the spatially variable detec- 
tor response, and mirror effective area (Bourdin et al. 2008). We 
estimated the ICM brightness from the spatial distribution of low 
energy X-ray photons (with energy lower than 2.5 keV), since 
their spectral distribution is expected to have a weak dependence 
on temperature. 

We used the X-ray data set in order to map the ICM temper- 
ature. To do so, we used the wavelet spectral-imaging algorithm 
proposed in Bourdin et al. (2004) and Bourdin et al. (2008). 
X-ray spectral-imaging can be performed only on data of suf- 
ficiently high count statistics, and it is necessary to gather pho- 
tons within large detector regions, the shape of which needs to 
be optimised to detect thermal structure. To overcome this major 
difficulty, we developed a strategy based on wavelet transforms, 
where we first estimated the ICM temperature with its associated 
confidence range at various angular scales and locations in the 
field of view. The detection of the temperature features and the 
reconstruction of an ICM temperature map was finally obtained 
by applying a threshold to the wavelet coefficients. The wavelet 
analysis was performed following a shift-invariant algorithm us- 
ing B -spline wavelet coefficients, and applying a threshold oiler 
confidence level above the noise fluctuation. Point sources have 
to be masked before building the map, to avoid pollution from 
the ICM signal at different spatial frequencies. 
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Table 1. XMM-Newton-EPIC observations used in our analysis, with effective exposure time corresponding to each instrument. In 
brackets: fraction of the useful exposure time after solar-flare "cleaning". 



Cluster name 


XMM-Newton 


Center coordinates 


MOS1 effective 


MOS2 effective 


PN effective 




obs. IDs 




exposure time (ks) 


exposure time (ks) 


exposure time (ks) 


A2933 


0305060101 


01h40m41.2s -54°33'26.0" 


21.5 (69.3 %) 


27.7 (89.0 %) 


16.3 (68.5 %) 


A2440 


0401920101 


22h23m52.6s -01°36'57.0" 


25.3 (55.4 %) 


24.9 (54.3 %) 


17.8 (50.3 %) 


A2384 


0101902701 


21h52ml4.2s-19°42'19.8" 


17.2 (65.4 %) 


16.8 (63.6 %) 


10.2 (51.6%) 



3. Optical observations: data description and 
reduction methods 

All three clusters were allocated time at ESO for wide-field 
imaging with the WFI instrument at the 2.2m telescope (to- 
tal of 9.75 hours), and extensive multi-object spectroscopy with 
EFOSC2 at the 3.6m telescope (total of 6 nights), in 2003, 2005, 
and 2007, respectively (programs 072.A-0595, 075.A-0264, and 
079.A-0425). As wide-field imaging of A2933 in service mode 
was not performed, we used the APM galaxy survey (Maddox 
et al. 1990) and the DSS images in the analysis of this cluster. 
We retrieved the relevant galaxy catalog in a region of 30'x30' 
covering A2933 from the Southern sky catalogue based on the 
UKST SES R survey. For A2384 and A2440, we obtained imag- 
ing of the central 30'x30' field in the R (filter ESO/844) and 
B (filter ESO/878) passbands. For each filter, 8 dithered images 
were obtained leading to a total exposure time of 40min. These 
sets of images were reduced and combined using the ESO/MVM 
package Alambic (Vandamme et al. 2002), and the catalogs were 
extracted with SExtractor (Bertin and Arnouts 1996). The typi- 
cal seeing is 1 arcsec. The magnitudes used in the present paper 
are the total magnitudes "MAG auto' provided by SExtractor. 
Stars and galaxies were classified (up to a magnitude of 21 
in both bands) by extracting the stellar sequences within the 
magnitude-half light radius diagrams. Fainter objects were all 
considered as galaxies. Finally, catalogs with B - R color were 
built by associating objects with positions in the B and R cata- 
logs differing by less than 1.5 arcsec. 

For A2440 and A2384, the final catalog of galaxies includ- 
ing both B and R passbands and covering a field of 30' x 30' was 
used to analyze the color properties and the spatial distribution 
of the cluster galaxies. The color-magnitude diagram was used to 
identify the red sequence defined by the cluster elliptical popula- 
tion (Lopez-Cruz et al. 2004). We selected galaxies with colors 
within 3 o~ of this relation as high probability members of the 
cluster and from this selection we derived the projected density 
maps for different magnitude cuts in the R-band (corresponding 
typically to L*+l, L*+2, L*+3). In the case of A2933, R band 
magnitudes from the APM catalog were used to compute density 
maps without color selection. 

Spectroscopy was performed with EFOSC2 using the grism 
#03, which covers the spectral range 3050-6010 A, with a 
FWHM of 7.5 A. In general, for each pointing we obtained 
2 x 45min exposures. Wavelength calibration was performed 
in real time, taking spectra of arc calibration lamps (Helium- 
Argon) after each exposure. Data reduction was performed using 
our dedicated IRAF package speXtra for automatic extraction 
and wavelength calibration of spectra. Radial velocities were ob- 
tained using the cross-correlation technique (Tonry and Davis 
1981) with the rvsao package. Cross-correlation was performed 
with velocity standards observed with the same instrumental 
configuration during the observing run. The typical velocity er- 
ror is ~ 50 km/s, an estimate that we could confirm by perform- 



ing a cross-check against the data available in the literature (see 
section 4.3). 

The spectroscopic catalogs resulting from our observations 
are listed for the three clusters in Tables |5j [6j and |7J where 
columns are as follows: 1) identification number; 2) and 3) right 
ascension and declination (J2000.0); 4) radial velocity; 5) veloc- 
ity error for the cross-correlation; 6) R TR parameter (Tonry and 
Davis 1981; when R TR > 3, the cross-correlation redshift can 
be considered as reliable); and 7) quality flag for the redshift (0: 
high precision, 1: medium precision). 

Objects with a velocity within a ±5000 km/s window cen- 
tered on the mean velocity were selected as cluster member can- 
didates. Possible interlopers were rejected by the gap technique. 
The velocity location and scale of the various clusters (and sub- 
clusters) were determined with the biweight estimator using the 
ROSTAT package (Beers et al. 1990). 

Ten normality tests (provided by ROSTAT) were applied to 
the data. The Dip test for unimodality (Hartigan and Hartigan 
1985) was computed with the diptest package implemented by 
Martin Maechler in the R environment. The P va iue was computed 
in each case according to the value of the Dip test and the num- 
ber of objects according to the table provided in Hartigan and 
Hartigan (1985). The results of this analysis are displayed in 
Table [2] where columns are as follows: 1) name of the cluster 
(or subcluster); 2) and 3): right ascension and declination of the 
brightest galaxy (BCG) of the (sub)cluster; 4) angular radius 6 
of the circle in which galaxies were selected for the ROSTAT 
analysis; 5) and 6) estimation of location and scale with the bi- 
weight technique; 7) number of redshifts used in the analysis; 
8) radial velocity of the BCG; 9) number of statistical tests ex- 
cluding Gaussianity at more than 10% confidence; 10) and 11) 
value of the Dip test and of its P- value; and 12) mean tempera- 
ture computed within radius 6 and its error. 

To test for multi-modality, we tried in each case to fit the 
entire velocity distribution with a mixture of Gaussian functions. 
For this purpose, we have used the powerful EMMIX algorithm 
(McLachlan and Krishnan, 1997; McLachlan et al.1999). This 
program is quite flexible and allows a variety of choices in the 
fitting. In particular, it does not require to introduce any guess for 
the initial partition of velocities and the form of the covariance 
matrix. We applied EMMIX to the data, trying to fit the velocity 
distribution with mixtures of N g Gaussian distributions, with N g 
varying from 1 to 5. We retained the best fit on the basis of the 
P-value, and in the case of identical values of P- value, on the 
criterion of the smallest partition number. When a mixture of 
gaussians had been fitted to the data, we recomputed the value 
of the location and scale of each partition with ROSTAT. 

4. A2933 

4.1. X-ray gas morphology and thermal structures 

In X-rays, A2933 is a well defined bimodal system. In Fig. [2] we 
show the temperature map of A2933 with the X-ray brightness 
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Fig. 1. Galaxy density maps as a function of galaxy luminosity. From top to bottom: A2933, A2440 and A2384. The field of view 
displayed is 15'x 15' for A2933 and A2440, and 30'x30' for A2384. From left to right, magnitude cuts are: R < 18, R < 19, R < 20 
and R < 21. The density maps of A2440 and A2384 have been computed for red sequence galaxies. 



contours superimposed. Looking at the low energy contours, 
both components look quite regular and have similar bright- 
nesses. We clearly detect the interaction zone as a hot region lo- 
cated in-between the two maxima. The northern unit (A2933N) 
is slightly elongated along the east-west direction, while the 
southern unit (A2933S) is orientated along the general NW- 
SE direction connecting the two units. In the temperature map, 
the most prominent feature is the hot region between the two 
units, which is statistically significant, because the exposure time 
of our XMM/EPIC observation was specifically adapted to en- 
sure good statistics in the region of interaction. The tempera- 
ture distribution of A2933N appears to be quite elongated, while 
A2933S looks more regular. 



4.2. Galaxy density distribution 

The morphology of Abell 2933 has been very poorly studied at 
optical wavelengths. In the following analysis, projected density 
maps were inferred in the field of the cluster from galaxy posi- 
tions and magnitudes in the APM catalog. 

The galaxy density distribution shows a bimodal structure, 
with two subclusters located on a NW/SE axis at a separation of 
~ 7 arcmin (Fig[T]). The relative importance of the two compo- 
nents varies with the magnitude cut-off. Selecting bright galax- 
ies (R < 18), the NW subcluster (A2933N) is dominant. When 
including fainter objects, the SE (A2933S) subcluster becomes 
progressively more prominent (Fig[T]). At faint magnitudes, the 
NW subcluster shows an extended tail in the NW direction. In 
Fig[3] we also show the isodensity contours corresponding to 



R < 20, superimposed on the DSS image of A2933. A2933N 
hosts a couple of very bright galaxies, and is centered on a BCG 
with a very close companion, while A2933S hosts a single bright 
elliptical galaxy. One can also note a small density clump at the 
southern extremity of A2933S. There is a general alignment of 
the two BCGs, the two subclusters, and the global cluster struc- 
ture along the NW/SE direction joining A2933N and A2933S. 



4.3. Galaxy velocity distribution 

A2933 was included in the Southern Abell Redshift Survey 
(SARS; Way et al. 2005), and on the basis of 53 redshifts Muriel 
et al. (2002) estimated the cluster global velocity and velocity 
dispersion (27 709 ±105 km/s and 759 ± 72 km/s, respectively); 
they also showed that the velocity dispersion as a function of ra- 
dius appears to be constant to 5 Mpc/h from the cluster center. 
Their data, however, sample the cluster and its environment on 
large scales and at relatively bright magnitudes. Since we are in- 
terested in the merging process, we observed the central region 
including the two main clumps at deeper magnitudes. 

We measured 71 redshifts in the 30' x 30' field of A2933, 
listed in Table [5] from which 47 objects are identified as clus- 
ter members. We have only 11 redshifts (among 53) in common 
with them, because of the different region sampled. Excluding 
one discrepant case, the redshifts of common objects are all in 
excellent agreement, with a mean velocity difference of 14 km/s 
and a standard deviation of 59 km/s. From our redshift sam- 
ple, we analyzed the velocity distribution in the central region 
of A2933. We measured a location of Cbi = 27 281 ± 103 km/s 
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Fig. 2. Left: Galaxy density map of A2933 (magnitude cut: R<20) overlaid on the X-ray luminosity contours (EPIC-XMM counts 
in the 0.5-2.5 keV band corrected for background and vignetting). The black crosses indicate the BCG positions. Right: ICM 
temperature map of A2933 overlaid on the X-ray luminosity contours (EPIC-XMM data analyzed by means of wavelet spectral- 
imaging; see Bourdin et al. 2004, 2008 for details). 



and a scale of S bi = 682 ± 65 km/s (Table [2|. Our estimate of 
velocity location is significantly lower (~ 400 km/s) than that 
found by Muriel et al. (2002). These differences may be caused 
by variations in sampling. Their sample covers a much larger 
area, that is not centered on the two subclusters but is shifted to 
the south, which we show to be a region of higher mean velocity. 
We checked that their estimate of mean velocity and our estimate 
for southern component (see below) are consistent. 

The velocity histogram (plotted in Fig [4]) clearly shows a 
multi-modal distribution, with a sharp peak at ~ 26 200 km/s, 
and a more dispersed structure, with two wider peaks at ~ 26 800 
km/s and ~ 27 500 km/s. In Fig [3] cluster members are indi- 
cated by different symbols and colors for galaxies with veloci- 
ties within each of the three peaks in the histogram. Most of the 
galaxies with velocities located in the sharp peak at ~ 26 200 
km/s (purple squares) are located in A2933N. Moreover, the re- 
gion of A2933S, including its southern extension, is populated 
mostly by galaxies in the highest velocity peak at ~ 27 500 km/s, 
while objects with velocities within the peak at ~ 26 800 km/s 
are mainly in the A2933N region. 

To quantitatively assess the significance of the multi- 
modality in the velocity distribution suggested by the histogram, 
we applied our set of statistical tests (running a set of tests is 
useful because they have different sensitivities to subclustering; 
see Pinkney et al. 1996). Ten out of 10 of the normality tests 
provided by ROSTAT do not find significant deviations from 
Gaussianity (see Table [2]). However, the Dip-test excludes uni- 
modality at more than 2% significance. 

We finally applied EMMIX, where the data were fitted 
without a priori constraints on the covariance matrix. A three 
Gaussian mixture inferred a P- value of 0.03, indicating that the 
Gaussian hypothesis is excluded at a very high significance level. 
This mixture fits the main structure with a Gaussian of mean 
V = 27 364 km/s and dispersion cr = 501 km/s (partition 1), 
the sharp peak at lower velocity with a Gaussian having mean 



V ~ 26 145 km/s and a very low dispersion cr = 65 km/s (par- 
tition 2), and the small excess at high velocities (partition 3). 
Results are listed in Table [3] The sharp low velocity excess con- 
sists of six galaxies: all but one are located in the northern sub- 
cluster (purple squares in Fig. [3]). These galaxies are likely to 
form a low velocity group probably infalling into A2933N. In 
the following, we exclude partition 3 (with only 2 objects) from 
the analysis. No significant results are obtained with EMMIX 
when fitting with mixture with a higher partitioning to the data. 



4.4. X-ray/optical combined analysis 

The projected density maps of A2933 for galaxies and for gas 
provide similar views of the cluster (Figj2j): the X-ray and op- 
tical isocontours of the NW and SE subclusters, A2933N and 
A2933S, closely correspond and are centered on the cold cores 
in the temperature map (Fig[2j. The X-ray contours are roughly 
centered on the BCGs for both the NW and SE subclusters. 
These results suggest that the system is in a pre-merger phase. 
However, in the case of the NW subcluster, the BCG seems 
slightly offset eastwards with respect to the X-ray centroid. The 
inner X-ray isocontours and the temperature map are elongated 
along a SE-NW axis connecting the two brightest galaxies of 
A2933N. 

As the density distribution of galaxies shows a clear bimodal 
behavior, we analyzed separately the velocity distributions of 
A2933N and A2933S, selecting the galaxies within a radius of 
3.5' from each subcluster center (see Fig[3]). The corresponding 
histograms are displayed in Fig|5] The velocity distributions of 
the two subclusters are offset from each other, the NW subclus- 
ter lying at lower velocities than the SE subcluster, but there is 
significant overlap. Moreover, they are quite different. A2933N 
has a multi-peak structure with the sharp excess at low velocities 
previously identified with EMMIX, while A2933S has a more 
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Fig. 3. DSS region (17'x20') centered on A2933 (North is 
up and East is to the left). Galaxies identified as cluster 
members from spectroscopy are marked with different sym- 
bols corresponding to the three major peaks in the veloc- 
ity histogram (purple squares: [26000,26500]km/s, cyan dia- 
monds: [26500,27200]km/s, orange circles: [27200,28500]km/s) 
which is plotted in Figj4j The isocontours corresponding to the 
galaxy density map with R < 19 are superimposed. The veloc- 
ity histograms for the galaxies in the two circles (which have a 
3.5 arcmin radius and are centered on the two subclusters) are 
displayed in Fig(5] 



continuous distribution with an extended low velocity tail that 
falls in the velocity range of A2933N. 

Departure from normality was tested when analyzing both 
subclusters individually, and neither the normality tests nor the 
Dip test allowed us to significantly exclude the hypothesis of a 
Gaussian velocity distribution. The ROSTAT analysis of the two 
velocity distributions shows that the scales (velocity dispersions) 
are comparable for both subclusters. These scales are in good 
agreement with those expected from the measured X-ray tem- 
peratures, assuming the typical scaling relation between <x and 
kT (Lubin and Bahcall 1993, Girardi et al. 1998, Wu et al. 1998, 
see Fig(6]), suggesting that gas and galaxies are in equilibrium. 

The dynamical analysis also confirms a significant offset for 
locations (mean velocities) of ~ 730 ±210 km/s: A2933S has 
a significantly higher mean velocity than A2933N (~ 27 648 
km/s vs. ~ 26919 km/s). This trend is also followed by the ve- 
locities of the brightest cluster members: the BCG in A2933N 




2.3xl0 4 2.4x10" 2.5x10* 2.6x10* 2.7x10* 2.8x10* 2.9x10* 3x10* 3.1x10* 
cz [km/s] 



Fig. 4. Velocity histogram of A2933 (binning of 250 km/s). The 
best Gaussian fit for the whole distribution (dotted line) is cen- 
tered on the vertical solid line which gives the location value. 
Location and scale of the Gaussian (i.e.mean velocity and ve- 
locity dispersion) were estimated with ROSTAT. We also show 
the two Gaussian functions corresponding to partition 1 and 2 
of the best mixture 3 partitions fit by EMMIX (dashed lines) 
and the composite function (solid line). 



.3x10* 




2.4x10* 2.5xl0 4 2.6 



<10* 2.7x10* 2.8xl0 4 2.9x10* 3x10* 3.1x10* 
cz [km/s] 



Fig. 5. Velocity histograms of A2933N (blue) and A2933S 
(red). For each subcluster we included galaxies within a ra- 
dius of 3.5' from its X-ray center. Solid lines: velocities of the 
brightest galaxy in the Northern (blue) and Southern (red) sub- 
clusters. 
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Fig. 6. The cr - T x relation for the subclusters A2933N 
and A2933S (red), A2440B and A2440C (green), A2384N 
and A2384S(blue). The straight line is the relation cr = 

10 2.47 ± 0.08 r 0.67 ± 0.09 (m et ^ 



(BCG1) and its close companion have comparable velocities 
(27 350 km/s and 27 424 km/s); the second brightest galaxy has 
a lower velocity (26951 km/s), while the brightest galaxy in 
A2933S (BCG2) has a higher velocity (28 028 km/s). This im- 
plies that there is a significant offset of ~ 600+ 140 km/s between 
BCG1 and BCG2. 

The velocities of the two BCGs also differ from the mean ve- 
locity of their host subcluster by 505 ± 165 km/s and 380 ± 200 
km/s for A2933N and A2933S, respectively, which is statisti- 
cally significant for A2933N. 

However, they coincide with the locations of the highest 
peaks in the velocity histograms of their respective subclusters 
A2933N and A2933S (excluding the sharp peak at ~ 26200 
km/s), and these peaks are slightly shifted from the mean ve- 
locity of each subcluster, due to the skewness of the velocity 
distributions (Fig [5]). All these results suggest that the velocity 
distributions have started to mix with each other, but that the 
structure of the two subclusters is still not strongly affected; this 
implies that the clusters have not yet crossed each other. 

Therefore the combined X-ray-optical analysis indicates 
that A2933N and A2933S are in an advanced pre-merger stage. 
The presence of two very bright galaxies along an NE-SW axis, 
the elongations of the inner X-ray contours and the temperature 
map along the E-W direction of A2933N, and the existence of 
a distinct low velocity, low dispersion velocity component, sug- 
gest that A2933N has also undergone some previous merging 
and is still accreting a small group. 



4.5. Dynamical analysis 

We now provide mass estimates for the subclusters and the to- 
tal mass of A2933 (analogous sections are devoted to A2440 
and A2384). For each system, we discuss the uncertainties in 
the mass estimates. An obvious uncertainty is related to the as- 
sumption of dynamical or hydrostatic equilibrium. For example, 
velocity dispersions might not simply reflect the subcluster po- 



tential and might be overestimated. Another uncertainty is due to 
the limited field in which it is possible to estimate the harmonic 
radius: the circle within which we estimate this radius cannot 
intersect the nearby subcluster(s). As a consequence, the virial 
radius might be underestimated. The mass may also be underes- 
timated due to the pressure term. 

4.5.1. Mass estimates 

From the previous analysis, A2933 is very likely in a pre-merger 
stage, in which subclusters begin to interact but cores have not 
yet crossed each other. We may then assume that the states of the 
two subclusters do not strongly deviate from dynamical equilib- 
rium, so that the virial theorem and the hydrostatic equilibrium 
can be safely applied to determine optical and X-ray mass esti- 
mates. These estimates will be used in the following to constrain 
the parameters of the collision by means of a two-body analysis. 

In the optical, we estimated the harmonic radius of the two 
subclusters from the projected distribution of all galaxies with 
R < 20 and within a radius of 3.5 arcmin from the BCG; from 
the harmonic radius we derived the virial radius and, in com- 
bination with the velocity dispersion, we estimated the virial 
mass. Finally, extrapolating the density profile, we estimated the 
mass M 2 oo corresponding to a density contrast p/p c = 200 (for 
more details see e.g. Maurogordato et al. 2008). We find for 
A2933N M 2 oo = 2.0 ± 0.4 x 10 14 M Q and for A2933S M 200 = 
(2.2 ± 0.4) x 10 14 M , with a ratio M A29 ^ S IM A29 ^ N = 1.1 ± 0.6. 

We can also estimate the mass ratio using the X-ray scaling 
relations. The mean X-ray temperatures of A2933N and A2933S 
are respectively kT A29 ^N = 2.1 1 ±0.25keV and £7^29335 = 2.7 ± 
0.2keV. Assuming M oc T^f 2 , we obtain a mass ratio between the 
two components around 1.5:1. 

On the other hand, when we estimate masses using the 
brightness and temperature profiles of each unit, assuming hy- 
drostatic equilibrium for each subcluster, we find that M^2933v = 
2.6 ± 0.4 x 10 14 M o and M A293 3S = 2.76 ± 1.5 x 10 14 M o leading 
to a mass ratio of §^ = 1.2 ± 0.6. 

MA2933N 

While the exact values could be systematically affected by 
deviations from equilibrium, there is a good consistency between 
the optical and X-ray mass estimates, and we can conclude that 
the two subclusters have comparable masses, with a total mass 
for A2933 of M 200 ~ 5 x 1O 14 M . 

4.5.2. Two-body analysis 

We applied the two-body dynamical formalism (Beers et al. 
1982; Gregory and Thomson 1984) to the two subclusters in 
A2933. One of the variants of this method establishes a rela- 
tionship between the total mass M of the bimodal cluster and the 
angle a between the plane of the sky and the collision axis of the 
two clumps (Barrena et al. 2002, Ferrari et al. 2005). The input 
parameters derived from the observations are the relative radial 
velocity between the two subclusters V r and their projected spa- 
tial separation R p . For these observables, we chose to take the 
values of the projected spatial separation and the radial relative 
velocity (at rest) between the two BCGs (R p = 0.75 Mpc and 
V r = 550 km/s), which are more clearly defined than that of the 
subclusters. We obtained similar results when using values from 
the subcluster centroids. In addition, the formalism requires us to 
define the time to that has elapsed since the epoch of the last en- 
counter between the two subclusters. We tested different scenar- 
ios, including a pre-merging scenario (to ~ 12.2 Gyr, i.e. the age 
of the Universe at the redshift of A2933, implying that the two 
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Fig. 7. The sum of the virial masses of A2933N and A2933S 
as a function of the projection angle a. The horizontal lines 
show the mass estimate (full line) and errors (dashed lines) 
(see Sect. 4.5. The projected distance of 0.75 Mpc and ra- 
dial velocity difference at rest of 550 km/s are derived from 
the two BCG's location and velocities. The two systems were 
at zero separation 12.2 Gyr ago. The dotted line represents the 
Newtonian criterion for gravitational binding. Two bound in- 
coming (Bla and Bib) and one bound outgoing (BO) are found 
compatible with our mass estimates. 



systems have never yet crossed each other), a recent (to = 0.5 
Gyr) and an older (to = 1.0 Gyr) post-merger event. Considering 
the total mass of the cluster derived from optical observations 
(5 x 10 14 Mq), the post-merger cases provide only one possi- 
ble bound outgoing solution (i.e. the two subclusters are moving 
apart). In the pre-merger case, Fig |7] reveals two possible bound 
incoming solutions, the first with a ~ 15° (Bla) and the second 
with a ~ 75° (Bib), and a bound outgoing solution (BO). For 
Bla, the "real" separation between the two subclusters appears 
to be small (~ 0.8 Mpc) and the relative velocity high (V ~ 2000 
km/s), while for Bib the two components appears to be more 
distant (~ 3.5 Mpc) and the relative velocity comparable to the 
projected observed one (V~ 560 km/s). The results are listed in 
TableH 

4.6. The scenario for A2933 

The combined X-ray optical analysis has shown: a) two subclus- 
ters, visible in both X-ray and optical maps, and associated with 
cool cores detected in the temperature maps; b) a hot region in 
between the two subclusters; c) a perturbed velocity distribu- 
tion, showing a velocity offset (at rest) of ~ 700 km/s between 
the two subclusters. This suggests that the interaction between 
the two subclusters is in its initial phase. This lead us to exclude 
the post-merger cases in the two-body analysis, and to examine 
the pre-merger one. Solutions with large values of a, Bib, and 
UO are quite improbable since they would imply large spatial 
separations between the subclusters (~ 3.5 Mpc and 11 Mpc, re- 
spectively), which is unlikely taken into account the clear signs 
of interaction between the two merging units. The most likely so- 
lution for A2933 is thus (Bla), i.e. a two-body pre-merger nearly 



on the plane of the sky (a ~ 15°). The high relative velocity 
(~ 2000 km/s) between the two components implied by this so- 
lution is a common feature in merging clusters. 



5. A2440 

Mohr et al. (1996) performed a combined X-ray and optical anal- 
ysis of A2440, based on optical imaging and spectroscopy (48 
redshifts), and Einstein X-ray data. Our work significantly ex- 
tends this study, being based on a larger spectroscopic sample 
(more than double the total number of redshifts) and new X-ray 
data, including the gas temperature maps. 



5.1. X-ray gas morphology and thermal structures 

The X-ray observations detect a gas emission with an elongated 
structure connecting two brightness peaks (B) and (C) to a north- 
ern and a less luminous component (A) (Fig. [8]). The two bright- 
ness peaks are associated with cool cores (kT ^ 2.8-3.5 keV), 
and are surrounded by hotter gas (kT ^ 5 keV); a significant es- 
timate of the temperature of the (A) component is not possible 
due to the strong contamination from the innermost cluster re- 
gions. 

The southern emission peak (C) is delimited by gas bright- 
ness and temperature jumps across the white sector shown in 
Fig[8] As discussed above, these discontinuities are signatures 
of a cold front delimiting the edge of the southern cold core. 

To investigate the nature of the gas density jump observable 
across the white sector in Fig [8] we model the gas 3D structure 
with disrupted density and temperature profiles (see Eqs. (20) 
and (21) of Bourdin et al. 2008). As shown in Fig[l0| fitting 
the slope, jump position, and amplitude of these functions re- 
veals a discontinuity in gas density, D n = 1.82 ± 0.06, and 3D 
temperature, Dj = 1.76^ ^, located at 110 kpc to the south of 
the cool core. The pressure continuity that is measurable across 
the jump (Dp = Dj/D n = 0.97 ± 0.05) identifies the discontinu- 
ity as a cold front. As already observed in various interacting 
systems, this cold front is likely to delimit a stripped cold core 
moving outwards from the hotter embedding ICM. 



5.2. Galaxy density distribution 

The galaxy density maps computed for RS galaxies for different 
magnitude cuts (Fig[T) show that the general structure of A2440 
is strongly elongated along a NE/SW axis and includes multi- 
ple clumps. At bright magnitudes (R < 18), four groups are de- 
tected: the two most significant subclusters B and C are identi- 
fied as the main cluster, and two fainter ones are detected at the 
NE extent: Al and A2 . When including faint galaxies, the two 
components corresponding to the main cluster B and C progres- 
sively merge to appear as a single dominant elongated structure 
(R < 21), while the NE groups Al and A2 are less visible and a 
new (but faint) component A3 is detected on the western side. 

The three brightest galaxies are located near the three density 
peaks Al, B, and C identified at R < 18 and R < 19 (Fig.pl). 

There is a strong alignment between the main axis of the 
bright galaxies, the position angle of the subclusters and the 
NE/SW axis of the whole cluster. 

5.3. Galaxy velocity distribution 

The first analyses of the dynamics of A2440 were performed 
by Beers et al. (1991) and Mohr et al. (1996). We now have a 
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Table 2. Characteristics of the clusters and subclusters. Location, scale, and X-ray mean temperatures of the different subclusters 
are computed in circles of 6 P arcmin centered on the position of the BCG (associated with the corresponding X-ray maximum). The 
measured velocity of the BCG is also listed. 
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Table 3. Results of the mixture of Gaussian with EMMIX. Errors have been estimated with ROSTAT. 
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Table 4. Two-body model solutions for the A2933N-A2933S, A2440A-A2440(B+C), and A2384N-A2384S systems. Various pos- 
sibilities are considered for the evolutionary phase: the pre-merger case (to = 12.2Gyr), and post-mergers seen at to after the first 
passage. For each solution, we derive the angle a between the line connecting the two components and the plane of the sky, the 
spatial separation of the subclusters R, their separation at maximum expansion R m and their relative velocity V. 



much larger dataset: from our EFOSC2 observations, 97 red- 
shifts have been obtained in the 30' by 30' field centered on 
A2440. We merged this catalog with redshifts available in the 
literature, which led to the final redshift catalog of 150 objects 
(with 10 twice observed objects). A comparison with 10 mea- 
sured redshifts available in the literature gives a mean difference 
of 30 km/s and a standard deviation of 60 km/s. Restricting this 
sample to cluster members, we are left with 103 objects. The 



velocity histogram is shown in Fig 12 



In Fig[TT] cluster members are indicated by purple squares 
and orange circles when belonging to one of the two major ve- 
locity peaks in the histogram at ~ 27 000 km/s and 28 000 km/s 
respectively. One can see that galaxies belonging to both veloc- 
ity peaks populate subclusters B and C, while all galaxies in sub- 
cluster Al belong to the highest velocity peak. 

Measurements of location and scale of the velocity distribu- 
tion with the biweight estimator give C B i = 27 251 + 93kms~ l 
and S B i = 940 + 7 Ohms- 1 . While none of the 10 ROSTAT nor- 
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Fig. 8. A2440: Left: Galaxy density map (mag. cuts: R<19) overlaid on the X-ray luminosity contours (EPIC-XMM counts in the 
.5-2.5 keV band corrected for background and vignetting). Black crosses indicate the BCG positions. Right: ICM temperature map 
overlaid on the X-ray luminosity contours (EPIC-XMM data analyzed through wavelet spectral-imaging, see Bourdin et al. 2004, 
2008 for details). 
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Fig. 9. Same as Fig.^for A2384. 
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mality tests identify significant deviations from Gaussianity, the 
Dip test rejects unimodality at the 1 percent level (Table [2]) . 

EMMIX fits a mixture with two partitions isolating the main 
component and an excess of two galaxies at V ~ 28 058 km/s 
with a very significant P-value (0.01). With three partitions, 
EMMIX isolates the main component with mean velocity V ~ 
27 080 km/s and velocity dispersion cr = 890 km/s (partition 1), 
a very sharp Gaussian (cr ~ 33 km/s), corresponding to the pre- 
viously mentioned excess with mean velocity V ~ 28 058 km/s 
(partition 2), and a group of two galaxies at ~ 30 337 km/s (par- 
tition 3), which are not considered in the following. The best- 



fit Gaussian functions (partition 1 and 2) are plotted in Fig. 12 
However, the relatively high value of P va iue obtained (0.23) in 
dicates that this three partition fit is less significant than the two 
partition one. No significant results are obtained with EMMIX 
when fitting with a larger number of partitions. 

5.4. A2440: X-ray/optical combined analysis 

Examining Fig [8] where we show the X-ray isocontours and the 
galaxy projected density distribution, we see that the central X- 
ray bimodal structure (B+C) is nearly coincident with the optical 
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Fig. 10. Top: Gas brightness and projected temperature profiles observable across the white sector in Fig. ^ for A2440. Bottom: 
ICM density and pressure profile estimated above the Southern gas clump in A2440. These profiles reveal a density discontinuity at 
constant pressure or "cold front" feature (see vertical dashed lines on the plots). 



bimodal structure (B+C) in the density map at R < 19 (and co- 
incident with that at R < 18). The X-ray A component is coinci- 
dent with the NE subcluster Al identified in the galaxy density 
maps, while the subclusters A2 and A3 are not detected in X- 
rays. The three X-ray maxima A, B, and C are well centered 
on the three corresponding BCGs (Fig [8]). We have seen that at 
fainter magnitudes the optical bimodal structure merges into one 
structure: its density peak is centered between the X-ray maxima 
B and C. These three subclusters, identified in X-rays and with 
optical counterparts, are referred to as A2440A, A2440B, and 
A2440C. 

To analyze the velocity distribution in A2440A, A2440B, 
and A2440C, we defined three subsamples including galaxies 
within circular regions of radius 2.5 arcmin (the largest compat- 
ible with negligible overlap) and centered on the corresponding 
X-ray maxima, which coincide with the BCGs positions. These 
circular regions are shown in FigfTT) 



In Fig 13 



we show the corresponding velocity histograms 
with different colors (blue, green, and red for A2440A, A2440B, 
and A2440C, respectively). The histograms indicate that most 
galaxies in A are in the velocity bin at ~ 28 000 km/s (EMMIX 
partition 2). Within the limits of the small number of counts per 
bin, there is no evidence of another significant segregation in 
velocity for galaxies in the regions B and C: the velocity distri- 
bution seems spatially mixed. 

The results of the ROSTAT analysis for the regions corre- 
sponding to the three X-ray subclusters are listed in Table [2] 
Both the mean velocity and velocity dispersions of the B and 
C regions are comparable, and consistent with the global cluster 
values. This suggests two possible alternatives: either A2440B 



and A2440C have not yet interacted, and the merging is occur- 
ring in the plane of the sky, or they have already crossed, but 
are seen just after core passage. The velocity of the BCG in 
the B region, BG2, is close to the velocity of its host subclus- 
ter, while BG3 in the C region has a (~ 2<x) velocity offset of 
(418 ± 188 kms~ l ). The velocity dispersions of the two main 
components A2440B and A2440C are higher than the values ex- 
pected from their X-ray temperature (see Fig(6]), but the devia- 
tion from the cr - T x relation is not statistically significant. 

The mean velocity of A2440A is « 150kms~ l higher than 
that of the main cluster at a 3cr level. The small velocity disper- 
sion of A is more typical of a group. The velocity of the brightest 
galaxy in A, VG1 = 27 925 kms -1 , is consistent with the mean 
velocity of its subcluster, implying that it is at rest in the poten- 
tial well. This dynamical analysis is in very good agreement with 
that obtained by Beers et al. (1991) and Mohr et al. (1996). 

These results show that A2440A is a group at higher veloc- 
ity with respect to the main system, which includes the two sub- 
clusters A2440B and A2440C: these two components must have 
recently crossed each other or are close to merging in the plane 
of the sky. 

We also applied the normality tests to the three subclusters. 
In the case of the A subcluster, the 10 ROSTAT tests and Dip 
test exclude normality at more than the 10% level. This result 
remains unchanged when excluding a low velocity galaxy that 
is 3000 km/s offset from the major velocity peak. We note how- 
ever that the analysis of this subcluster relies on only 11 (10) 
velocities. For subclusters B and C, all the normality tests (ex- 
cept for one in the case of B) and the Dip test are consistent with 
a unimodal Gaussian distribution. 



12 



S.Maurogordato et al.: Merging history of three bimodal clusters 




15 | < < i < | < < < < | < < < < | < < < < 



Fig. 11. WFI R-band image of the A2440 field (12'xl4'). 
Galaxies identified as cluster members from spectroscopy (Flag 
and 1) are marked with different symbols according to the ve- 
locity range (purple squares: [24000,27600] km/s; orange circles 
[27600,30000] km/s). The isocontours of red sequence galaxy 
density maps with magnitude limit R < 19 are superimposed. 
The velocity histograms of the three regions with 2.5 arcmin ra- 
dius delimited by circles (corresponding to the three subclusters 
A2440A, A2440B, and A2440C) are displayed in Fig[l3] North 
is up and East is to the left. 



5.5. Dynamical analysis 
5.5.1. Mass estimates 

We estimated the mass of the three components in A2440 follow- 
ing the same method as for A2933. However, to determine the 
harmonic radius we selected the galaxies in the red sequence, de- 
fined using our B and R photometric catalogues, with R < 20; in 
this way, background contamination was minimized. For each 
component, we selected galaxies satisfying the above criteria 
within a radius of 3.0 arcmin from the respective BCG. 

As a caveat, we emphasize that we consider subclusters po- 
tentially in a post-merger phase, in which case one expects im- 
portant distortions in morphology and velocity field. However, 
previous analyses have shown that the velocity distributions in 
both components B and C can be assumed to be Gaussian. We 
therefore assume that the two main components, A2440B and 
A2440C, are now not far from dynamical equilibrium. In con- 
trast, A2440A has a non-Gaussian velocity distribution, but this 
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Fig. 12. Velocity histogram of A2440 (binning of 250 
km/s). The best Gaussian fit for the whole distribution (dot- 
ted line) is centered on the vertical solid line which gives 
the location value. We also show the two Gaussian func- 
tions (dashed lines) corresponding to partitions 1 and 2 in 
the best 3 partition mixture model by EMMIX and the com- 
posite function (solid line). 
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Fig. 13. Velocity histograms of A2440A (blue), A2440B 
(green) and A2440C (red). For each subcluster we in- 
cluded galaxies within a radius of 3' from its center. Dotted 
lines: velocities of the brightest galaxies in A2440A (blue), 
A2440B (green) and A2440C (red). 
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was determined using only a few redshifts. However, its charac- 
teristics are similar to those of a group, and its contribution to 
the whole mass should be a priori negligible. 

We find for A2440A, A2440B, and A2440C, respectively, 
that M 2 oo = 0.1 ± 0.02 x 10 14 M , M 200 = 5.6 ± 1.1 x 10 14 M , 
and M200 = 5.4 ± 1.1 x 10 14 M©. As previously discussed, B and 
C appear to be the main subclusters; the estimated total mass of 
A2440 is M A244 o - 1.2 X 10 15 M Q . 

5.5.2. A2440: Two-body model 

The two-body model was applied to the central system com- 
posed of A2440B and A2440C, with input parameters R p = 
0.442 Mpc and V r = 200 km/s derived from the BCGs. As 
expected in terms of both projected separation and velocities, 
the two-body model is not well constrained by the observations. 
Bound solutions exist in both the pre-merger case with to = 12.2 
Gyr (two incoming and one outgoing), and in the post-merger 
case (one outgoing for to < 0.5 Gyr, two incoming, and one out- 
going for to > 0.5 Gyr). Therefore, according to the two-body 
analysis, A2440B and A2440C may be systems seen either be- 
fore or after first encounter. 

We also tested the two-body system composed of A2440A 
and the whole complex A2440(B+C). To estimate R p and V r , 
we used the position and velocity of the centroid of the whole 
system A2440(B+C), obtaining Rp = 0.6 Mpc and V r = 730 
km/s. In the pre-merger hypothesis (to - 12.2 Gyr), there are 
two bound incoming (Bla, Bib), and one bound outgoing (UO) 
solutions (Fig. 14). If seen after first passage, one bound out- 
going solution exists (BO). The different solutions are listed in 
TableH 



5.6. The scenario for A2440 

We can now use the complete set of information collected from 
our analysis of A2440 to determine the merging scenario. The 
bimodal X-ray emissivity of the main component is evidence of 
the interaction of two massive subclusters being likely to un- 
dergo or to have undergone a merging event. At first sight, the 
existence of two units with cool cores, and a close correspon- 
dence between the gas and bright galaxy distributions, which are 
well centered on the BCGs, are indicative of a pre-merger event. 
This was also the conclusion of Mohr et al. (1996) for A2440, 
who observed the similarity between galaxy projected density 
maps and X-ray emission maps from Einstein. 

However the pre-merger hypothesis is strongly disfavored by 
the lack of any strong temperature enhancement between the two 
maxima, as we would expect in a pre-merger event (RS01). The 
similar velocity distributions of the two subclusters can be un- 
derstood in terms of both a pre-merger and post-merger, but with 
different implications: in the pre-merger case it is indicative of 
a merging in the plane of the sky, while in the post-merger case 
it is due to a very recent collision. The observed difference of 
clustering for faint and bright galaxies, with a good coincidence 
between the gas and bright galaxies centroids for A2440-B and 
A2440-C subclusters, but a segregation between the centroids of 
the gas and faint galaxies distributions also indicate that merging 
has already occurred. These results have already been found for 
several post-merger clusters (Biviano et al. 1996; Maurogordato 
et al. 2008). 

The evidence of a cold front delimiting the core of the 
Southern subcluster also favors a post-merger scenario, where 
the stripped subcluster is now moving away from its compan- 




Fig. 14. Two-body model for the system A2440A/A2440(B+C). 
The horizontal lines show the mass estimate (full line) and errors 
(dashed lines) (see section 5.5. 1| ). The projected distance of 0.6 
Mpc and the radial velocity difference at rest of 730 km/s were 
derived from the position and velocity of the two centroids. The 
two systems are assumed to have been at zero separation 12.2 
Gyr ago. The dotted line represents the Newtonian criterion for 
gravitational binding. Two bound incoming (Bla and Bib) and 
one unbound outgoing (UO) solutions are found compatible with 
the mass estimates. 



ion, as shown for instance in A1201 by Owers et al. (2009). 
The velocity offset of BG2 with respect to A2440-C also implies 
that there is some dynamical activity in this region. To make this 
scenario compatible with the observed bimodal structure, we ar- 
gue that the collision must have occurred with a relatively high 
impact parameter because a too small impact parameter would 
have led to the disruption of the two former subcluster cores 
during the violent relaxation of the gas. The similar values ob- 
tained for the temperatures and velocity dispersions of A2440B 
and A2440C suggests that these components have comparable 
masses. If one identifies the hot structures surrounding A2440B 
and A2440C as remnants of a shock wave, a scenario with an 
equal mass collision with a high impact parameter after the max- 
imum core collapse, when the shock wave begins to propagate 
towards the outskirts of the new structure, is therefore very likely 
to explain the A2440(B+C) complex. 

The A2440A subcluster is probably a group, with a signifi- 
cant velocity offset (~ 750 km/s) with respect to the main clus- 
ter A2440(B+C). The coincidence of the BCG position with the 
maximum of the gas and galaxy distribution and the consistency 
of the BCG velocity with the subcluster mean velocity, suggest 
that the BCG is at the center of the subcluster in dynamical equi- 
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librium. We should expect a Gaussian velocity distribution for 
this subcluster, which is not the case: this may be due to the 
small number of available redshifts (10) and contamination from 
galaxies not belonging to A2440A. 

The lack of signatures of gas compression between A2440A 
and A2440B excludes a pre-merger hypothesis implying small 
physical separations (Bla). The two viable pre-merger solutions 
are then incoming (Bib) or outgoing (BO), both with a major 
component along the line of sight (a ~ 80 degrees). The post- 
merger case remains unlikely, as there is a good correspondence 
between the gas and galaxy distributions and the X-ray temper- 
ature map does not show any evidence of shocked gas in the 
periphery of A2440A. Therefore A2440A is probably a group 
infalling for the first time into the main cluster component along 
the NE/SW axis. 

6. A2384 

6.1. X-ray gas morphology and thermal structures 

Looking at Fig|9| it is clear that the X-ray morphology of A2384 
is very peculiar. The general shape of the whole cluster is very 
elongated, including a northern primary maximum, connected 
to a southern secondary maximum by a continuous gas distri- 
bution. The temperature of this gas is around 3-4 keV, while its 
surroundings are hotter (up to 4-5 keV). Both subclusters exhibit 
cool cores. Another interesting feature in the temperature map is 
the 5 keV region spreading along the eastern region. 



The hypothesis of a massive system is improbable, given the low 
cluster optical luminosity and richness. The battery of normality 
tests do not indicate significant departures from normality. The 
Dip test excludes unimodality with a P- value of 0.1. 

cluster members in the three velocity peaks are 



In Fig 15 



identified with different symbols and colors. We have very few 
redshifts in A2384S, because it is a poor structure. A striking 
feature is the high fraction of objects in the filament joining 
A2384N and A2384S, which are associated to the second ve- 
locity peak (~ 28 500 km/s). 

For A2384, the best-fit solution was found by EMMIX for a 
partition of three Gaussians roughly centered on the peaks vis- 
ible in the velocity histogram, with a P- value of 0.16. This in- 
dicates that, while the fit quality is greatly enhanced by using 3 
Gaussians instead of 1, the null hypothesis of Gaussianity is not 
rejected. The three Gaussian functions are visualized in Fig[l6j 
the sum of these functions reproduces the data quite well. The 
results are summarized in Table [3] The whole velocity distri- 
bution can be reproduced by the combination of two Gaussians 
with intermediate velocity distributions (~ 365 km/s and ~ 665 
km/s corresponding to partitions 1 and 3, respectively) at mean 
velocities ~ 1200 km/s higher and lower than the mean velocity 
of the whole cluster, and by a third strongly peaked component 
of low dispersion, cr ~ 250 km/s peaked at ~ 28 500 km/s (cor- 
responding to partition 2). In the following section, we discuss 
the spatial distribution of the galaxies assigned to the different 
velocity partitions with respect to the X-ray and optical density 
distribution. 



6.2. Galaxy density distribution 

A2384 is one of the best examples of colinear distribution of 
substructures in a cluster (see West, Jones, Forman 1995). At 
all magnitude limits, the galaxy distribution shows a very elon- 
gated structure, extending over 1.2 Mpc along the N-S axis 
(Fig [I]); embedded in this structure, there are two main sub- 
clusters roughly centered on two bright galaxies (which we call 
BCG1 and BCG2 in the following; see Fig 15). The northern 
subcluster is more densely populated than the southern one. At 
bright magnitudes (R<18), the galaxy distribution is strongly 
concentrated in the northern (A2384N) component. Besides 
BCG1, in A2384N there are two other very bright galaxies, 
while in A2384S there are no bright galaxies around BCG2. 
Near A2384N, a clump of galaxies to the west causes an elonga- 
tion in the density map. While A2384S is very elongated along 
the general N/S axis of the system, the internal contours of 
A2384N are oriented towards the NW/SE, and rotate to the W-E 
and general NE/SW direction on larger scales. When including 
fainter objects, the isodensity contours of the two components 
are elongated towards each other (see Fig[T]). 

6.3. Galaxy velocity distribution 

Before our observations, very few redshifts were available in the 
literature for this cluster. In our 30' x 30' field, only 4 galax- 
ies already had a redshift; we reobserved two of them, and our 
redshifts are in agreement with the literature values within the 
estimated errors. Our final catalog includes 84 redshifts, and the 
number of cluster members is 56. 

The mean location of the cluster, CBI= 28263 ±154 km/s, is 
well defined, but the scale is quite high (SBI= 1 1 14 ± 120 km/s). 
The velocity distribution (Fig 16) is very broad and has three 
main peaks. This extended velocity distribution could be either 
due to a very deep potential (massive system) or multi-modality. 



6.4. A2384: X-ray/optical combined analysis 

Looking at Fig|9j both the gas and the galaxy maps approx- 
imately define the same bimodal structure. The northern and 
southern overdensities (A2384N and A2384S) correspond to the 
cool cores in the X-ray temperature map. The X-ray maximum 
of A2384S is coincident with the BCG2 position, while BCG1 
is slightly offset westwards (15 arcsec) with respect to the X-ray 
maximum of A2384N (see Fig|9]). This offset cannot be due to 
astrometric errors, as the precision of the optical images is su- 
perior to 1", while the precision of X-ray images is superior to 
4", and there is a very good correspondence between the optical 
and X-ray positions of point sources. As BCG1 also has a sig- 
nificant velocity offset with respect to A2384N (see below), this 
indicates that it is not at rest in the subcluster potential. However, 
the velocity distribution of A2384N has a large dispersion, and 
at the cluster redshift 15 arcsec corresponds to a physical pro- 
jected separation of about 25 kpc, which means that BCG1 is 
close to the cluster center, which is probably included in its ex- 
tended halo. 

As already noted, the velocity dispersion of A2384 is unex- 
pectedly high with respect to its optical luminosity and to the low 
values of its X-ray luminosity and temperature. A cluster with 
velocity dispersion ~ 1000 km/s should have an X-ray tempera- 
ture T ~ 7 KeV, according to the cr - T x relation (Fig(6]), which 
is much higher than the observed value. Moreover, from a weak 
lensing analysis Cypriano et al. (2004) fit a velocity dispersion 
of 737 ±126 km/s or 797 ±108 km/s depending on the assumed 
profile (SIS and SIE, respectively). These velocity dispersions 
corresponds to T ~ 3.45 keV and T ~ 4.02 keV, respectively, 
which are consistent with the observed values in X-rays. This 
suggests that the measured velocity dispersion is overestimated. 

We analyzed the velocity distribution in A2384N and 
A2384S, selecting galaxies in circles of radius 5 arcmin cen- 
tered on the X-ray maxima. The velocity histograms in the cor- 
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Fig. 15. WFI R-band image (20'x25') centered on A2384. Galaxies identified as cluster members from spectroscopy (flags and 
1) are marked with different symbols as in the previous cases (purple squares: [24000,28000] km/s, cyan diamonds: [28000, 29000] 
km/s and orange circles: [29000,31000] km/s). Isocontours corresponding to red sequence galaxy density maps at R < 19 are 
superimposed. The velocity histograms of the three regions delimited by circles with 3 arcmin radius, which include the two 
subclusters and the intermediary region, are displayed in Fig [18] North is up and East is to the left. 



responding regions are shown in Fig 17 The mean velocity of 
A2384S is ~ 675 km/s higher than that of A2384N. The same 
trend is followed by the BCGs of the two subclusters, which 
show an even larger offset of ~ 1000 km/s. We also found that 
the position of the BCG1 has a small offset with respect to the 
centroid of A2384N, while the position of BCG2 corresponds 
to the centroid of A2384S. The velocity of BCG1 is 600 km/s 
lower than the mean velocity of galaxies in A2384N, while the 
velocity of BCG2 is consistent with that of A2384S (see Table 

The velocity dispersions of the two subclusters are also quite 
large: ~ 1200 km/s for A2384N and ~ 900 km/s for A2384S. 
Both subclusters are above the cr - T x relation; the deviation is 
at the 3<x level for A2384N, but within the errors in the case of 
A2384S. ROSTAT normality tests reject the Gaussian hypoth- 
esis for neither A2384-N, nor A2384S. The Dip-test however 
excludes unimodality for A2384-N (A2384-S) with a P- value of 
0.1 (0.01). This suggests that the subclusters are still dynami- 
cally perturbed. 

We also defined three subsamples by selecting galaxies 
within a radius of 3 arcmin from the centers of A2384N, 
A2384S, and the intermediate region (see Fig [18]). The galaxies 
in the northern component show a large spread in their veloc- 



ity distribution, which is skewed towards lower values and has a 
peak at ~ 28 500 km/s, while galaxies in the southern and espe- 
cially in the intermediate region are more strongly concentrated 
in the same peak; this peak corresponds to EMMIX partition 2. 



6.5. Dynamical analysis 



6.5.1. Mass estimates 



We applied a method similar to that applied to A2440, determin- 
ing the harmonic radius of each component by selecting only the 
galaxies in the red sequence, with R < 20 and within a radius of 
5 arcmin from the respective BCG. 

However, as we previously discussed, the velocity distribu- 
tion of this system is not simple. Our mass estimates could there- 
fore be severely affected. We find for A2384N and A2384S that 
M 2 oo = 1.4 ± 0.3 x 10 15 M and M 200 = 0.6 ± 0.1 x 10 15 M , re- 
spectively. As the velocity dispersion is probably overestimated, 
we consider in the following M 2 oo = 2.0 x 10 15 M as an upper 
bound to the sum of the mass of the two components. 
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Fig. 16. Velocity histogram of A2384 (binning of 250 km/s). 
The best Gaussian fit for the whole distribution (dotted line) 
is centered on the vertical solid line which gives the location 
value. Location and scale of the Gaussian were estimated with 
ROSTAT. We also show the three Gaussian functions (dashed 
lines) corresponding to the best mixture model by EMMIX and 
the composite function (solid line). 



6.5.2. Two-body model 

We applied the two-body model to A2384N and A2384S, assum- 
ing as values of R p and V r the projected separation and relative 
velocity at rest of the BCGs (Rp = l.UMpc, Vr = 1000 km/s). 
We tested several scenarios: the pre-merger case (to = 12.2 Gyr), 
and a range of post-mergers (to = 0.2, 0.5, 1., 1.5 and 2 Gyr). In 
the pre-merger case, two incoming bound solution and one un- 
bound solutions were found. In the post-merger case, only an 
unbound solution exists if the merging is recent (to < 0.5Gyr), 
and a bound outgoing solution exists if the merger is older 
(to > 0.5Gyr). Figure 19 shows the case of a post-merger seen 
1.0 Gyr after the encounter for which a unique bound outgoing 
solution exists. At to > l.5Gyr, two bound incoming solutions 
and a bound outgoing solution are possible. 



6.5.3. A2384 proposed scenario 

The bimodal structure of A2384 detected in both the galaxy and 
gas density maps suggests that the two components either are 
about to or have undergone a merger event. The presence of two 
well-separated components with cold cores, coincident in the X- 
ray and the optical, is often observed in the case of pre-mergers. 

However, the trail of gas and the bridge of galaxies extending 
from A2384N to A2384S are difficult to explain if the subclus- 



Fig.17. Velocity histograms of A2384N (blue), and A2384S 
(red). For each subcluster we included galaxies within a radius 
of 5' from its center. Velocity bins of 500 km/s are used. Dashed 
lines: velocities of the brightest galaxies in A2384N (blue) and 
A2384S (red). 



ters have not yet interacted. If A2384 were a pre-merger, its tem- 
perature map should exhibit a compression region between the 
two subclusters, which is not the case. The velocity distribution 
shows a large dispersion (cr ~ 1200 km/s), inconsistent with the 
value of the temperature if it simply reflects the mass of the sys- 
tem. This high value of the velocity dispersion and the velocity 
distributions of the two components being mixed imply that the 
merging has already occurred. A2384N shows several signatures 
of previous dynamical activity: a slight segregation between gas 
and galaxies, an offset of the position of the BCG from to the 
peak of the X-ray emitting gas, and of its velocity from the mean 
velocity of the subcluster. This suggests that the dynamical state 
of A2384N is still disturbed possibly due to a previous collision 
with A2384S or/and a secondary merger event within A2384N 
itself. 

At first sight, A2384S is more regular than A2384N. It has 
a good spatial coincidence between its gas and galaxy distribu- 
tions, coincidence between the position and velocity of the BCG 
and the X-ray center and mean velocity of the subcluster, respec- 
tively. Its velocity distribution is very broad and non-Gaussian, 
but this may be due to contamination by objects in the periphery 
ofA2384N. 

This complex structure could be understood in a scenario in 
which a low mass cluster of galaxies has crossed the environ- 
ment of a more massive one and has been stripped of a large 
fraction of its gas and galaxies. The "head" of the system would 
then correspond to A2384S, and the tail of galaxies belonging to 
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Fig. 18. Velocity histograms of A2384N (blue), the intermedi- 
ate region (green) and A2384S (red). For each subcluster we in- 
cluded galaxies within a radius of 3' from its center. Velocity 
bins of 500 km/s are used. Dashed lines: velocities of the bright- 
est galaxies in A2384N (blue) and A2384S (red). 



the original cluster affected by A2384N would lie in the inter- 
mediary region. 

We conclude that the most likely scenario for A2384 is a 
post-merger between two unequal mass clusters. However, if the 
merging has already occurred, we have to explain why the two 
cool cores have remained unaffected by the merging process. 
Numerical simulations (Poole et al. 2006) have shown that the 
initial cool cores of the primary and secondary components sur- 
vive the first core crossing, and disappear after the second peri- 
centric passage. 

It is however difficult to constrain the parameters of the sce- 
nario, as the mass of A2384 is not well determined. From the 2- 
body analysis, bound solutions exists only if the merger event is 
older than 0.5 Gyr. Older mergers (1-1.5 Gyr) with larger phys- 
ical separations between the units are favoured by the lack of 
signatures of compression in the X-ray gas (Table [3]). For merg- 
ers older than 1.5 Gyr, one can see that incoming solutions, cor- 
responding to second infall, are again possible. However, these 
also predict small values of the real separation between the two 
components, which are unlikely from the lack of signature of 
compression of the gas in the temperature maps, as noted be- 
fore. One is then left with the most probable hypothesis of two 
bound outgoing subcluster s seen more than 1.0 Gyr after the first 
passage, with a collision near the line of sight. Figure 17 shows 
as an example the bound outgoing solution in a 1.0 Gyr post- 
merger, which is in good agreement with both our X-ray and 
optical observations. 



Fig. 19. The sum of the virial masses of A2384N and A2384S as 
a function of the projection angle a. The horizontal lines show 
the mas s val ue (full line) and errors (dashed lines) as estimated 
in Sect. |6.4| The projected distance of 1.14 Mpc and the radial 



velocity difference at rest of 1000 km/s were derived from the 
two BCGs positions and velocities. It is assumed that the two 
systems are 1.0 Gyr after the first passage. The dotted line rep- 
resents the Newtonian criterion for gravitational binding. One 
unbound outgoing (UO) is found compatible with our mass esti- 
mates. 



7. "Idealized" simulation & XMM emulation 

In the previous sections, we combined optical and X-ray obser- 
vations to constrain the merging scenario for each one of our 
three bimodal systems. To help define the scenario and test its 
consistency, we tried to reproduce the observed X-ray proper- 
ties with numerical simulations. Adopting an approach similar 
to RS01, we placed two clusters in a box and allowed them to 
evolve under gravity in an adiabatic framework. The initial pa- 
rameters of the collision are inferred from the X-ray and optical 
observational constraints. 

To this aim, we developed a set of programs that allow us to 
reproduce the observational process and analyze the simulated 
clusters in the same way as we analyze real ones, adopting an 
approach similar to e.g. Gardini et al. (2004) for XMM observa- 
tions. 

Using an emission plasma model, where gas pressure and 
density are known everywhere, we can compute the emitted X- 
ray photons as a function of telescope area and exposure time. 
Fixing the direction of observation, we collect photons from the 
different cells along the line of sight, taking into account that we 



18 



S.Maurogordato et al.: Merging history of three bimodal clusters 



observe clusters in projection (as long as we can assume that the 
thin plasma approach is valid). 

The first step in developing a good "idealized" merger sim- 
ulation is to check the way we implement the hydrostatic equi- 
librium in each (sub)cluster (hereafter "unit") to be sure that, 
during collision, what we see is really due to the collision and 
not a spurious result of numerical errors in the initial units. Once 
this is achieved (by allowing a unit to evolve alone in the middle 
of the box for a very long time), we are able to build collision 
simulations. 

The total mass profile of our "perfect cluster" follows a clas- 
sic NFW profile, Pdmi?) oc r( - l+ l /rs ^ 2 (Navarro et al. 1996), while 
the gas follows a profile described in Suto et al. (1998) adapted 
to keep the gas mass fraction constant at large radii. The dark 
matter profile is derived from the difference between the total 
mass profile and the gas. We note that the NFW profile is de- 
rived from pure dark matter simulations, thus this approach is 
likely to be safe as long as the gas mass fraction is not too high 
(thus not in the central region of cooling-flow clusters). 

The velocity dispersion of dark matter particles is com- 
puted (as in RS01) using the virial equation -^:lpdm(r)cr 2 (r)] = 

GM{r) , v 
—PdmV)> 

After fixing the gas-to-mass ratio, the concentration param- 
eter C200 = ' m e mass ratio between units and the impact pa- 
rameter (defined as the minimal distance between the two initial 
trajectories in unit of the NFW scaling radius r s ), we reconstruct 
a pre-collision scenario by computing two units as described 
above. We place these two units in a rather large box to avoid 
any border effects. We compute the initial velocity considering 
the pre-merger evolution as a free fall of the two units, exactly 
in the way described in RS01. We then compute the evolution of 
the system with the AMR hydro-NBody simulation RAMSES 
(Teyssier et al. 2002). At each time step, for all the cells in the 
simulated volume, we store the pressure, velocity, and density of 
the gas, and position and velocity of dark matter particles. After 
running these simulations, the following steps are the computa- 
tion of the photons emitted by the simulated volume and their 
ray-tracing using the EPIC/XMM instrumental response. These 
steps are fully described in Bourdin et al. (2004). We finally have 
a dataset that can be reduced in the same way as real XMM/EPIC 
observations. 

All these previously described steps are important for a cor- 
rect comparison between observation and simulation. Because 
of the complexity of this procedure, it is impossible to explore 
the whole parameter space when determining the optimal col- 
lision parameters. By looking simultaneously at the brightness, 
temperature maps of the X-ray gas, and galaxy density maps, 
we first guess a set of possible parameters that we place into the 
complete simulation stream and iterate manually to optimize the 
set of parameters. Therefore, we do not pretend to fit the data, but 
we attempt to obtain a scenario reproducing the main features of 
the observed merging system. 

8. Comparison with simulations and discussion 

Taking as input the scenario inferred from observations for the 
main merger event, we obtained density and temperature maps 
as they would be observed with XMM, and compared them with 
the real observations. 

In the case of A2933, the most likely scenario is a pre-merger 
of two components of comparable mass for which there is a 
small angle between the collision axis and the plane of the sky. 
We simulated an equal-mass merger in the plane of the sky with 
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Fig. 20. Simulated wavelet reconstructed temperature map with 
luminosity contours superimposed, derived from the simulation 
which best reproduces the observed features of A2933. 



an impact parameter set to zero, and allowed it to evolve with 
time. The resulting temperature map and luminosity contours 
for this scenario are presented in Fig 20 The X-ray properties, in 
particular the temperature contrast between the high temperature 
region and the mean temperature of each unit and its relative an- 
gular distance, are closely reproduced by a pre-merger scenario 
~ 200Mvr before the core collapse. 

For A2440, which is expected to be a post-merger case, we 
simulated a variety of equal mass mergers with different impact 
parameters and tried to find the best configuration reproducing 
simultaneously the observed temperature and luminosity distri- 
butions. In Fig|2T] we present our best candidate. This scenario 
corresponds to a collision with an impact parameter of 15 r s seen 
450 Myr after the maximum core collapse. We managed to re- 
produce the most striking features revealed by the observations, 
in particular the hot region surrounding the X-ray maximum of 
each unit, and the temperature contrast between the hottest re- 
gion and the colder ones. Since our code follows a purely adi- 
abatic approach, we do not expect of course to reconstruct the 
two cool cores. 

In the case of A2384, we performed different trials allow- 
ing the impact parameter value to vary. If the impact parameter 
is not zero, it should be large enough to allow the gas of the 
smaller unit to follow the influence of its original potential well 
after the collision (we note that this is also necessary to explain 
the remaining "cool" core of each unit). At the same time, the 
impact parameter should be small enough to enable the gas to 
be stripped efficiently. In these limits, the gas is stripped along a 
rather significant distance as observed. Simulations by Poole et 
al. (2006) show a gas surface density similar to that of A2384, in 
the case of a 3:1 mass ratio seen at the first apocentric passage 
of the secondary cluster, at ~ 1.5 Gyr after the closest approach 
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Fig. 21. Simulated wavelet reconstructed temperature map with 
luminosity contours superimposed, derived from the simulation 
which best reproduces the observed features of A2440. 



(see their Fig. 4 central panel). In these scenarios, it is neverthe- 
less very difficult to keep the axial symmetry observed in lumi- 
nosity and in temperature maps and impossible to explain that 
the BCG of the south unit still lies at the X-ray maximum. 

On the other hand, in a null impact parameter collision, if a 
small group encounters a heavier unit, its gas can escape signif- 
icantly after the first maximum core collapse only if its central 
density is very high compare to the central density of the mas- 
sive object. Simulations by Poole et al. (2006) are able to strip 
gas even with a null parameter thanks to their initial gas pro- 
file (S oc r 11 ) profile, which leads to a high inner gas density 
in the core. In contrast, if the small unit has a concentration pa- 
rameter higher than the massive one (say by a factor 2) then the 
gas appears sticked to its own unit even after the first maximum 
core collapse. Unfortunately, we cannot simulate such a complex 
scenario in the adiabatic scheme since we have no cooling thus 
the inner temperature of the small unit increases strongly with 
the concentration parameter, preventing the correct description 
of the gas luminosity and temperature during the collision. 

9. Conclusion 

We have performed a combined X-ray/optical analysis of three 
bimodal clusters at low redshift, selected as merger candidates 
from the sample of Kolokotronis et al. (2001). We have con- 
firmed the disturbed dynamical state of these clusters and re- 
constructed the general trends of their main merging scenario. 
A2933 has been shown to most likely be an equal-mass ad- 
vanced pre-merger (~ 200Myr before the core collapse). A2440 
is a recent equal-mass merger (~ 450 Myr after core collapse) 
with another subcluster infalling along the main axis of the clus- 
ter, and A2384 is the result of an older collision between two 
units of mass ratio 1:3 (~ 1.0 - 1.5 Gyr after core collapse). 



We emphasize the complementarity of the X-ray and optical 
data in determining the parameters of the scenario (angle, mass- 
ratios, epoch) that were refined by running simulations of ide- 
alized cluster collisions. Our analysis confirms the efficiency of 
the selection based on the comparison of gas and density maps to 
identify merging clusters when (as is frequently the case) infor- 
mation on temperature and velocity distribution is not available. 
Moreover, the level of segregation between galaxies and gas pro- 
vides some indication of the merging stage. This approach is far 
more powerful than the density maps taken individually in ei- 
ther X-rays or the optical. From X-ray density maps only, the 
three bimodal clusters presented in this paper would have been 
expected to be in the very last phase before the first maximum 
core collapse, while two of them are shown to be post-mergers 
with surviving cores in the two units. The properties of the tem- 
perature maps are essential to assess the merging stage: the hot 
region in between the two sub-clusters in A2933 is indicative of 
a pre-merger, the spectacular hot structures in the periphery of 
A2440 and A2384 have been identified as remnants of a post- 
merger shock wave, and the cold front detected in the S/W re- 
gion of A2440 is consistent with a a post-merger. On the other 
hand, redshift information is needed to understand the dynamics 
of the system, determine the mass of the subclusters and the rel- 
ative motion of the BCGs, and finally date the merger event and 
constrain its axis using the two-body analysis. 

Numerical simulations have allowed us to refine the colli- 
sion parameters. Thanks to our XMM simulation tools, we cor- 
rectly handle all projection and instrumental effects. Thus we can 
rely on these "observed simulations" and select the most likely 
scenario based on measured criteria such as distances, tempera- 
ture, and luminosity. We have been able to reproduce the main 
features of A2933 and A2440, simulating the collision of two 
systems initially in complete hydrostatic equilibrium (with an 
approach similar to RS01) and adding a complete set of pro- 
grams to reproduce X-ray observations for simulations and ob- 
tain EPIC/XMM synthetic observations (with an approach simi- 
lar to Gardini et al. 2004). The case of A2384 is too complex to 
be described with our approach, but we have been able to place 
constraints on its collision parameters. 

The results obtained from the analysis of these clusters can 
be combined to our previous work on other MUSIC clusters to 
draw a more general picture of merging clusters properties. In 
the three clusters studied in the present paper, we have shown 
the existence of preferential axes following the general posi- 
tion angle of the cluster, as previously shown in other merging 
clusters (Arnaud et al. 2000; Plionis et al. 2003; Maurogordato 
et al. 2006). Subclusters generally host one or more BCGs. In 
most cases, the different merging subclusters and their bright- 
est BCGs are aligned along this direction. This corroborates the 
property of preferential alignment of BCGs with their host clus- 
ter (Niederste-Ostholt et al. 2010) on the scale of subclusters. In 
A2384, we have detected a spectacular filament of galaxies and 
gas probably stripped from the colliding group along the merg- 
ing axis. 

Another interesting aspect is the position and motion of the 
BCGs relative to their subclusters. We have found that their 
angular coordinates generally correspond (within the errors) to 
the X-ray centroids of their host subclusters (A2933S, A2440A, 
A2440B, A2440C, A2384S) but in some cases there is a small 
offset (A2933N, A2384N). Similarly, in some cases the radial 
velocities of the BCGs are consistent with the mean velocity of 
their respective subcluster (A2440A, A2440B, A2384S), but in 
other cases they are offset by 300-500 km/s (A2933N, A2933S, 
A2440C, A2384N). These offsets are signatures of ongoing dy- 
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namical activity in the subclusters, as found in A3921B (Ferrari 
et al. 2005) and in A2163A (Maurogordato et al. 2008), which 
were both identified as recent mergers. In an analogous way, 
the beginning of interaction between A2933N and A2933S, the 
probable recent mergers occuring within A2933N and A2384N, 
and the recent merger between A2440B and A2440C can explain 
the observed BCG offsets. 

We have also found that, in addition to the main merging 
event, a large fraction of MUSIC clusters appear to contain sec- 
ondary merging events (A2933N and A2384N) and infalling 
groups,e.g., A2440 (this work), A521 (Ferrari et al. 2003), and 
A2163 (Maurogordato et al. 2008). 

These results are consistent with a hierarchical scenario of 
structure formation where clusters form by successive mergers 
and accretion of matter along large-scale filaments, which may 
cause alignments with structures on various scales (Basilakos et 
al. 2006; Lee and Evrard 2007; Faltenbacher et al. 2008). 

Moreover, we have also shown that in recent post-merger 
clusters, such as A2440 (this work) and A2163 (Maurogordato et 
al. 2008), galaxies exhibit a strong luminosity segregation, sim- 
ilar to the case of the Coma cluster studied by Biviano et al. 
(1996). We have also found a deviation from the cr - T x rela- 
tion for the subclusters in a post-merger stage, the largest one 
(at the 3cr level) being detected for A2384N. At variance, both 
pre-merging subclusters A2933N and A2933S follow the cr-T x 
relation. 

The impact of the merging process on galaxy properties in 
these clusters, in particular on star formation, will be addressed 
in a forthcoming paper. 

To test the validity of these properties on a large sample, we 
are currently extending this work to a subsample of the C4 SDSS 
cluster sample (Miller et al. 2004) with both optical and X-ray 
available data. 
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Table 5. Radial velocity measurements in the field of Abell 2933. The full catalogue is available in electronic form at www.edpsciences.org 



Galaxy 


Right Ascension(J2000) 


Declination (J2000) 


v (km s l ) 


6 (km s l ) 


Rtr 


Flag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


1 40 14.73 


- 54 32 36.2 


184 


19 


16.8 





2 


1 40 16.47 


-54 29 51.7 


26159 


104 


2.7 





3 


1 40 17.54 


-54 28 15.3 


27500 


121 


3.3 





4 


1 40 18.80 


- 54 28 45.6 


27303 


80 


4.3 





5 


1 40 20.24 


-54 29 13.4 


26356 


101 


3.8 





6 


1 40 22.09 


- 54 29 28.6 


152 


50 


5.6 





7 


1 40 22.15 


- 54 29 28.6 


137 


21 


16.0 





8 


1 40 22.98 


- 54 32 10.6 


26623 


33 


10.5 





9 


1 40 25.52 


- 54 28 54.4 


26489 


66 


4.2 





10 


1 40 25.84 


_ 54 30 49.4 


26018 


106 


3.0 
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Table 6. Radial velocity measurements in the field of Abell 2440. The full catalogue is available in electronic form at www.edpsciences.org 



Galaxy 


Right Ascension(J2000) 


Declination (J2000) 


v (km s l ) 


6 (km s l ) 


Rtr 


Flag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


22 23 21.33 


- 1 40 54.0 


58154 


54 


3.5 


1 


2 


22 23 22.01 


- 1 42 31.2 


-136 


13 


18.4 





3 


22 23 23.56 


- 1 42 19.5 


27293 


27 


7.2 





4 


22 23 24.82 


- 1 41 38.1 


-134 


57 


4.7 


1 


5 


22 23 27.64 


- 1 40 57.3 


26389 


14 


0.0 





6 


22 23 28.80 


- 1 42 28.6 


-51 


11 


23.0 





7 


22 23 30.06 


- 1 41 50.0 


26670 


19 


15.1 





8 


22 23 30.75 


- 1 40 44.5 


92248 


41 


5.6 





9 


22 23 31.92 


- 1 42 21.2 


27415 


22 


13.4 





10 


22 23 33.70 


- 1 39 58.8 


26784 


54 


5.2 
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Table 7. Radial velocity measurements in the field of Abell 2384. The full catalogue is available in electronic form at www.edpsciences.org 



Galaxy 


Right Ascension(J2000) 


Declination (J2000) 


v (km s l ) 


6 (km s l ) 


Rtr 


Flag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


21 52 2.86 


- 19 41 36.4 


26810 


35 


7.6 





2 


21 52 3.13 


- 19 36 57.7 


18 


41 


8.6 





3 


21 52 4.20 


- 19 39 20.8 


-69 


23 


16.3 





4 


21 52 4.45 


- 19 30 17.7 


27250 


120 


0.0 


1 


5 


21 52 5.01 


- 19 32 45.3 


476 


94 


3.5 





6 


21 52 5.44 


- 19 36 24.9 


29289 


57 


6.2 





7 


21 52 6.10 


- 19 41 33.0 


27472 


25 


10.0 





8 


21 52 6.36 


- 19 35 54.6 


27454 


67 


4.5 





9 


21 52 6.71 


- 19 39 25.7 


28590 


88 


3.3 





10 


21 52 6.72 


- 19 42 43.1 


28927 


85 


2.8 


1 



